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In event by event hydrodynamic model, we have simulated 30-40% Au+Au collisions at RHIC 
and computed the slope parameter from the invariant pion distribution. In each event, the slope 
parameter fluctuates azimuthally. Fourier expansion coefficients T n for the slope parameter and 
the Fourier expansion coefficients v n for the azimuthal distribution 4^ are found to be strongly 
correlated. Strong correlation between the two expansion coefficients suggests that in addition to 
azimuthal distribution, fluctuations in the slope parameter of the invariant distribution can as well 
be used to study the final state momentum anisotropy in relativistic energy heavy ion collisions. If 
measured experimentally, they can serve as additional constraint for hydrodynamical modeling. 

PACS numbers: 47.75.+f, 25.75.-q, 25.75.Ld 


I. INTRODUCTION 

In recent years, there is much interest in event-by¬ 
event hydrodynamics. Event-by-event hydrodynamics 
takes into account that in nucleus-nucleus collisions, par¬ 
ticipant nucleon positions fluctuates from event to event. 
Effects of such fluctuations are most prominent on the 
azimuthal distribution of the produced particles. In a 
non-zero impact parameter collision between two iden¬ 
tical nuclei, the collision zone is asymmetric. Multiple 
collisions transform the initial asymmetry into momen¬ 
tum anisotropy. Momentum anisotropy is best studied 
by decomposing it in a Fourier series, 
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(j) is the azimuthal angle of the detected particle and ijj n is 
the plane of the symmetry of initial collision zone. In cen¬ 
tral rapidity region, for smooth initial matter distribution 
(obtained from geometric overlap of density distributions 
of the colliding nuclei), plane of symmetry of the colli¬ 
sion zone coincides with the reaction plane (the plane 
containing the impact parameter and the beam axis), 
ip n = \k j Rp,V?r. The odd Fourier coefficients are zero by 
symmetry. However, fluctuations in the positions of the 
participating nucleons can lead to non-smooth density 
distribution, which will fluctuate on event-by-event basis. 
The participating nucleons then determine the symmetry 
plane ( ippp ), which fluctuate around the reaction plane 
[l[. As a result odd harmonics v±, V 3 , v$ ..., which in cen¬ 
tral rapidity region are exactly zero for smoothed initial 
distribution, can be developed. In RHIC and LHC en¬ 
ergy collisions, odd harmonics has indeed been measured 
in experiments 0 [1] 4] 0 ■ Explicit event-by-event hydro- 
dynamic simulations of relativistic heavy ion collisions 
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While momentum anisotropy is best studied by Fourier 
expansion of the azimuthal distribution, in the present 
paper we explore the possibility of using another ob¬ 
servable, namely the slope parameter of the transverse 
momentum distribution for the produced particles. In 
experiments azimuthal angle integrated transverse mo¬ 
mentum distribution ( , d N , —) of produced particles are 
routinely measured. Angle integrated momentum distri¬ 
bution is of exponential nature. Slope (T) of the distribu¬ 
tion can be interpreted approximately as the temperature 
of the fireball produced in the collisions. In event-by¬ 
event hydrodynamics the initial medium is granular and 
even in central collisions have azimuthal dependence. As 
a consequence the freeze-out surface can have azimuthal 
dependence. The invariant particle spectra will depend 
on the azimuthal angle. Hydrodynamic simulations in¬ 
dicate that the invariant spectra (— d , N ,, ) remains 

v \dyp T dp T d(p J 

exponential even as a function of the azimuth. Asymme¬ 
try in the momentum distribution then be reflected as 
the azimuth dependent slope parameter {T(cj 1 )). Exper¬ 
imental measurements can be easily extended to study 
of azimuthal fluctuations in the slope parameter. The 
data can serve as additional constraint for hydrodynam¬ 
ical modeling. 

In the present paper, in event-by-event hydrodynam¬ 
ics, we have simulated ^^^=200 GeV 30-40% Au+Au 
collisions and studied azimuthal fluctuations of the slope 
parameter. When averaged over many events, azimuthal 
fluctuation in the slope parameter is very similar to that 
of azimuthal distribution; though less pronounced. In 
a single event however, the slope parameter appears to 
fluctuates more strongly over the azimuth than the az¬ 
imuthal distribution. In analogy to the Fourier expan¬ 
sion of the azimuthal distribution 4A. we have Fourier 

dtp ’ 

expanded azimuth dependent slope parameter. The ex¬ 
pansion coefficients of the slope parameter and that of 
azimuthal distribution appear to be strongly correlated. 
Strong correlations between two expansions coefficients 
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suggest that like the azimuthal distributions, azimuth 
dependent slope parameter can be used to characterise 
the momentum anisotropy. 

The paper is organised as follows: in section [XT] the 
hydrodynamic model used presently is briefly discussed. 
In section m simulation results are discussed and finally 
in section nv] conclusions are drawn. 

II. HYDRODYNAMIC MODEL 


tained by assuming that each participant deposits energy 
£o in the transverse plane, 

Npart 

e(x,y) m £ 0 ^2 S(x - x i: y - yi) (4) 

i= 1 

For use in fluid dynamical models, the discrete density 
distribution is smoothed by smearing the participant po¬ 
sitions by a Gaussian function, 


In event-by-event hydrodynamics, one generally solves 
for the energy-momentum and baryon number conserva¬ 
tion equations, 


dfjT^ = 0 , ( 2 ) 

= 0 (3) 

where T^ = (e-| -p)u IJ ’u l '—Pg ,J ’ 1 ' is the energy-momentum 
tensor and J M = ubu m is the particle 4-current, e, p, 
ns and u are energy density, pressure, net baryon den¬ 
sity and hydrodynamic 4-velocity respectively, = 

diag( 1,-1,—1,-1) is the metric tensor. Presently we 
assume that fluid is baryon free and disregard Eq(3] 
We also disregard any dissipative effect. Assuming 

boost invariance we solve the equations with the code 
AZHYDRO-KOLKATA pj] in (r = y/t 2 - z 2 , x, y, r/ s = 
\ In jz 2 -) coordinate system. Hydrodynamics equations 
(EqU]) are closed with an equation of state (EoS) p = 
p(e). Presently, we use an equation of state where the 
Wuppertal-Budapest mm lattice simulations for the 
deconfined phase is smoothly joined at T = T c = 174 
MeV, with hadronic resonance gas EoS comprising of all 
the resonances below mass m res = 2.5 GeV. Details of the 
EoS can be found in jUJ. 

Solution of hydrodynamic equations requires to specify 
the fluid energy density distribution £(x,y), velocity dis¬ 
tribution v x {x,y),v v (x,y) at the initial time. A freeze- 
out prescription is also needed to convert the informa¬ 
tion about fluid energy density and velocity to invari¬ 
ant particle distribution. We assume that at the initial 
time Tj=0.6 fm initial fluid velocity is zero, v x (x,y) = 
v v (x,y) = 0. The freeze-out temperature is fixed at 
Tp=130 MeV. For the initial energy density distribu¬ 
tion we resort to Monte-Carlo Glauber model. Details 
of the Monte-Carlo Glauber model can be found in !22j . 
In a Monte-Carlo Glauber model, according to the den¬ 
sity distribution of the colliding nuclei, two nucleons are 
randomly chosen. They are assumed to interact if the 
transverse separation is below For RHIC energy 

collisions <jnn ~ 42m&, a value used in the present sim¬ 
ulations. Transverse position of the participating nucle¬ 
ons (which are known in each event) will fluctuate from 
event-to-event. If a particular event has N part partici¬ 
pants, participants positions in the transverse plane can 
be labeled as, (aq, J/i), {x 2 , y 2 )--(x Npart , 2/jv part )- Energy 
density distribution in the particular event can be ob- 
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e(x,y )=£ 0 ^ gGauss{x - y,Xi,y - yi,cr) (5) 
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In f23|[24|], influence of the Gaussian width i.e. the 
smoothing parameter ex, on flow coefficients were stud¬ 
ied in Au+Au collisions. Elliptic and triangular flows 
are minimally influenced by the smoothing parameter a. 
Higher flow coefficients, however are influenced by the 
choices of a. Presently, we have used <r=0.25 fm. 


III. RESULTS 

A. Feasibility of slope parameter as a probe for 
momentum asymmetry 

Let us first demonstrate the feasibility of slope parame¬ 
ter in characterising momentum anisotropy of final state 
particles. We have simulated 30-40% Au+Au collision 
at RHIC energy, with two initial condition IC-1 and IC- 
2. IC-1 is a randomly chosen MC Glauber model initial 
condition and IC-2 is the initial condition obtained by 
averaging over average of 500 MC events. In FigQJi and 
(b), energy density distributions for the two initial con¬ 
ditions are shown. When averaged over large number of 
events, the initial condition is rather smooth as it is in 
smooth hydrodynamics (obtained from optical Glauber 
model calculations). Initial condition IC-1, correspond¬ 
ing to a single MC event clearly shows granular struc¬ 
tures. 

In FigJTJi* and d, angular dependence of 7r~ spectra 
from the evolution of initial distributions IC-1 and IC-2 
are shown. Spectra are shown for azimuthal angle </>=0, 
1.26, 2.51, 3.77, 5.03 and 6.28 respectively. In the inset 
of'CEfc and d, spectra are shown in a limited pt range. For 
IC-2, as expected, the spectra do not show any angular 
dependence. The angular dependence however is mani¬ 
festly presented in IC-1. We also note that at any angle 
the spectrum is of exponential nature, exp(—rriT/T). 
Azimuth dependent slope parameter then can be used 
to probe momentum anisotropy of the final state parti¬ 
cles. In smooth hydrodynamics, the azimuth indepen¬ 
dent slope parameter can be interpreted as effective tem¬ 
perature of a single fireball emitting particles. Azimuth 
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FIG. 1: (color online) (a) Distribution of initial energy density 
in a single MC event (IC-1), (b) distribution of energy den¬ 
sity averaged over 500 MC events (IC-2). (c) Transverse mo¬ 
mentum distribution for 7r“ from hydrodynamic simulation 
of initial condition IC-1, (d) Transverse momentum distribu¬ 
tion for 7r _ from hydrodynamic simulation of initial condition 
IC-2. 



0 1 2 3 4 5 6 



FIG. 2: (color online) (a) Variation of the fluctuations in the 
slope parameter AT((f)) with azimuth. The red line is for a 
single MC Glauber event and the filled circles with error bars 
are for average of 500 MC Glauber events, (b) same as in (a) 
but for fluctuations in the azimuthal distribution AN. 


dependent slope parameter then can be interpreted as 
the temperatures of multiple sources. 


B. Azimuthal variation of slope parameter in 
Event-by-event hydrodynamics 

For 500 MC Glauber events we have simulated ir~ pro¬ 
duction in 30-40% Au+Au collisions at RHIC. For each 
event we find the slope parameter ( T{q !>)) between pr= 1- 
2 GeV, and compute the fluctuations in T(cj>) as, 

T(d>) -T — IT dd>T(d>) 

A TM = . T = Jo £ W (7) 


Figj2jr, the red line is the fluctuations in the event IC-1. 
T{cj )) more oscillates more rapidly than the event aver¬ 
aged value. Fig0o shows the azimuthal variation A N(4>). 
Event averaged (ATV) fluctuations are similar to that of 
the slope parameter, though more pronounced. In Fig0o, 
the red line indicate azimuthal variation of AN in the 
event IC-1. In a single event, unlike AT, AN oscillates 
less rapidly. 


C. Fourier expansion of slope parameter 

As in EqJTj azimuthal anisotropy of slope parameter 
can be studied by Fourier expansion, 


In simulations we also have calculated the azimuthal 
distribution dA. j n each of the event. Fluctuations in 
azimuthal distribution are calculated in a similar way, 

dN _ ~dW 

AN (6) = d<t> _ ; 

dN_ 
d(j> 

In Fig 0 azimuthal fluctuations AT(</>) and AN((f>) are 
shown. In Fig0i, the circles with error bars shows the 
azimuthal variation for the event averaged fluctuations 
(AT(</>)}. When averaged over all the events, the slope 
parameter does not show rapid variation with azimuth. 
However, large error bars do indicate that at any az¬ 
imuth, event-by-event variations are rather large. In 


dN fTd^% 
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T(4>) = T 0 


1 + 2 ^2 T n cos{ruj) — nip n ) 

n 


,n = 1,2,3... 


(9) 

where the expansion coefficients are labeled as Ti, T 2 .. 
etc. In Eq]9] (also in EqQ]) ij} n is the participant plane 
angle. In event-by-event hydrodynamics, one generally 
characterise the asymmetry of the initial collision zone 
in terms of various moments of the eccentricity ( e n ) 


e n e 


inijjn 


f f e(x,y)r n e m ^dxdy 
II e(x, y)r n dxdy 


n= 1,2,3..(10) 


which also determine the participant plane angle 1 jj n . 
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FIG. 3: (Color online) Pictorial depiction of particle distri¬ 
bution with finite vi and vz in the transverse plane. The red 
lines in the left and right panels show the polar plot of the 
azimuthal distribution 4A. with finite V2 and finite V3. The 
black circles shows the distribution when V 2 =V 3 = 0 . 

To understand the implications of the coefficients T n , 
let us remember that the physical implications of the flow 
coefficients v n in Eq[l] In Fig[3]schematic diagram of the 
azimuthal distribution with finite V 2 and with finite 
U 3 are shown in polar coordinates. For finite V 2 , in the 
transverse plane, the distribution is elliptical, particles 
are preferentially produced at 0 = 0 and 0 = tt. The 
same for distribution with finite V 3 however is triangular 
in shape. Preferential emission occurs in three direction, 
0=0, 2tt/3 and 47 t/3 . In Fig[3] solid circles shows the az- 
imuthally symmetric distribution when V 2 or V 3 are zero. 
Microscopically, due to multiple collisions, initial asym¬ 
metry in spatial density distribution gets converted into 
asymmetry in the momentum space distribution. If it is 
assumed that particles are emitted from a hot fireball in 
the momentum space, then for finite V 2 (^ 3 ), momentum 
space density distribution of the fireball, in the trans¬ 
verse plane is elliptical (triangular). The coefficients v n 
in EqQ]are called flow coefficients as they indicate prefer¬ 
ential flow of particles in certain directions. Physical im¬ 
plication of the coefficients T n in Fourier expansion of the 
slope parameter is similar to that of the flow coefficients 
v n . In a distribution with finite T 2 , in the transverse 
plane fireball temperature is not uniform rather elliptical 
in shape, more along 0=0 and n and less in other direc¬ 
tions. Similarly fireball temperature will be of triangular 
shape in distribution with finite T 3 . In a fireball picture 
such angular dependence implies that a number of fire¬ 
balls participate in particle production. Considering the 
similarity with flow coefficients v n , in the following T n 
will be called temperature flow coefficients. Similar to 
flow coefficients, microscopic origin of temperature flow 
coefficients is also the asymmetry in the initial (spatial) 
density distribution. 

In a hydrodynamic model flow coefficients v n are re¬ 
sponse of the initial (spatial) asymmetry of the medium 
produced in the collisions. They are expected to be 
strongly correlated with initial eccentricity measured e n . 
Indeed, in smooth hydrodynamics elliptic flow is strongly 
correlated with initial eccentricity £ 2 - In event-by-event 


hydrodynamics correlation between flow coefficient 
with initial eccentricity measures e n has been studied pre¬ 
viously 0 , 0 > 0 , [29( 1 . It was shown that in event-by¬ 
event hydrodynamics, elliptic flow is strongly correlated 
with initial eccentricity. Comparatively weak correlation 
was observed between triangular flow and initial triangu¬ 
larity. Correlation between higher flow coefficients and 
asymmetry measures gets even weaker. De-correlation of 
higher order flows could be understood as due to nonlin¬ 
ear mixing of modes [l 6 [ . For example, Gardim et al jT|| 
showed that in order to correctly predict V 4 and V 5 , one 
must take into account nonlinear terms proportional e 2 
and £263 respectively. 


TABLE I: Correlation measure for (i) flow harmonics (v„) 
and eccentricity (e„), (ii) temperature flow coefficients (T„) 
and eccentricity (e n ), and (iii) between temperature flow co¬ 
efficients (T n ) and flow coefficient (u n ). 


'measure 



n— 1 

n = 2 

n = 3 

n = 4 

n = 5 

{y n ■> €n) 

0.656 

0.974 

0.846 

0.682 

0.715 

{Tn, e„) 

0.572 

0.797 

0.796 

0.651 

0.651 

{Tn , Vn) 

0.897 

0.824 

0.940 

0.926 

0.909 


As we have argued here, temperature flow coefficients 
T n are also responses of the initial (spatial) asymmetry of 
the system. How they are correlated with initial eccen¬ 
tricity measures is an interesting study. In Fig[|J we have 
studied the correlation between initial eccentricity (e„) 
and temperature flow coefficients (T„). For each of the 
500 MC-Glauber events the slope parameter T(0) was 
Fourier expanded to find out the temperature flow coeffi¬ 
cients T n , n=l,2,..5. In the left panel of Figj4l for the 500 
events, the temperature flow coefficients T n (n=l,2,..5) 
are plotted against initial eccentricity measure e„. For a 
perfect correlation T n oc e ra and all the events should lie 
on a straight line. It is evident from our simulations that 
temperature flow coefficients are weakly correlated with 
initial eccentricity measures. 

The result is not a complete surprise. From our previ¬ 
ous studies, we know that with the exception of elliptic 
flow V 2 , the flow coefficient v n are also weakly correlated 
with initial eccentricity measures. It is evident from the 
middle panel of FigjT] where we have plotted the flow 
coefficients v n against the eccentricity measures e n . As 
expected, only V 2 appears to be strongly correlated with 
£ 2 . Other flow coefficients are weakly correlated with 
eccentricity measures. In the right panel of FigQJ we 
have plotted temperature flow coefficients T n against the 
flow coefficient v n . Interestingly, temperature flow coef¬ 
ficients T n and flow coefficient v n appear to be strongly 
correlated. 

in [umi a quantitative measure was defined to quan¬ 
tify the correlation between flow coefficients and initial 
spatial asymmetry measure. 
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30-40% collision 
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FIG. 4: Each dot in the plots represent event-by-event hy¬ 
drodynamic simulation for flow coefficients v n , temperature 
flow coefficients T n . Left panel: Correlation between temper¬ 
ature flow coefficients T n with initial eccentricity measure e n . 
Middle panel: Correlation between flow coefficients v n with 
initial eccentricity measure e„. Right panel: Correlation be¬ 
tween temperature flow coefficients T n with flow coefficients 

Vn- 


c, 


measure 


(n) = ! - 


V n ,st.line(£n)] 

YsiHandom ( e ) “ ^Ume(<0] 2 


( 11 ) 


Cmeasure essentially measures the dispersion of the 
simulated flow coefficients from the best fitted straight 
line, relative to completely random flow coefficients. It 
varies between 0 and 1. If flow coefficients are perfectly 
correlated then v n oc e n and C me asure is identically unity. 
For completely random flow coefficients, C rneaS ure= 0- To 
obtain an even ground for comparison of Cmeasure for dif¬ 
ferent flow coefficients, the flow coefficients (v n ) and the 
asymmetry parameters (e n ) are scaled to vary between 0 
and 1. A similar equation can be used to quantify the 
correlation between temperature flow coefficients T n and 
initial asymmetry measure e n or between temperature 
flow coefficients T n and flow coefficients v n . In tabled] 


we have listed the Cmeasure values for correlation between 
v n and e n . between T n and e n and between T n and v n . 
Quantitatively compared to flow coefficients v n , tempera¬ 
ture flow coefficients T n are less correlated with e n . The 
difference is most prominent for second harmonic n=2 
and lessened in higher order harmonics. Cmeasure values 
also indicate that the correlation between temperature 
flow coefficients T n and flow coefficient v n are rather high 
for all n. T n and v n are much better correlated than the 
correlation between v n and e n or between T n and e n . 

Event-by-event hydrodynamic simulation result that 
the temperature flow coefficients T n and flow coefficients 
v n are rather strongly correlated strongly suggests that 
apart from the azimuthal distribution , the azimuthal 
fluctuations in the slope parameter of invariant spec¬ 
tra can as well be used to study final state momentum 
anisotropy. The fluctuations can be measured experi¬ 
mentally and can serve as additional constrains for hy¬ 
drodynamic modeling. 


IV. SUMMARY AND CONCLUSIONS 


In event-by-event hydrodynamical model for heavy ion 
collisions, due to fluctuations in participant positions, 
the initial conditions fluctuates event-by-event. These 
fluctuations are manifested as momentum anisotropy in 
the final particle distribution and in general are stud¬ 
ied by Fourier expanding the azimuthal distribution, 
the expansion coefficients characterising the momentum 
anisotropy. By explicit event-by-event hydrodynamic 
simulation we have shown that in event-by-event hydro¬ 
dynamics, the slope parameter of the invariant distribu¬ 
tion will also fluctuates with azimuth and can as well 
be used to study the momentum anisotropy. For 500 
MC Glauber events, we have simulated Au+Au collisions 
at RHIC and from the invariant pion spectra computed 
the azimuthal distribution ^ and from an exponential 
fit to the spectra between pt= 1-2 GeV, also the slope 
parameter T(tp). ^ and T(</>) were Fourier expanded 
and studied its correlation between expansion coefficients 
and initial eccentricity measures. With the exception of 
n=2, expansion coefficients v n for the azimuthal distribu¬ 
tion ^ are weakly correlated with the initial eccentric¬ 
ity measures e n - Without any exception, the expansion 
coefficients T n for the slope parameter are also weakly 
correlated with the initial eccentricity measures e n . The 
coefficients v n and T n are however strongly correlated. 
Strong correlation between v n and T n suggest that they 
complement each other. If measured experimentally fluc¬ 
tuations in slope parameter can be used to study momen¬ 
tum anisotropy and can serve as additional constraint for 
hydrodynamical modeling. 
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